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We discuss the B-meson light-cone wavefunction relevant for QCD factorization approach for 
exclusive Z?-meson decays. We derive the operator product expansion for the Z?-meson light-cone 
wavefunction, taking into account the local composite operators of dimension less than 6 and cal- 
culating the radiative corrections at order for the corresponding Wilson coefficients. The result 
embodies peculiar UV and IR behaviors of the Z?-meson light-cone wavefunction, the Sudakov- 
type double logarithmic effects and the mixing of the multiparticle states with additional gluons 
inside the B meson. The former effects are induced from the cusp singularity in the radiative 
corrections, while the latter is manifested by the participation of the higher-dimensional operators 
associated with the nonperturbative structure of the B meson. 
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For the exclusive B-meson decays such as B ^ nn, py,-.-, systematic methods have been 
developed using QCD factorization based on the heavy-quark limit [|l|]. Among the building blocks 
in the corresponding factorization formula of the decay amplitude, essential roles are played by 
the light-cone wavefunctions (LCWFs) for the participating mesons, which include nonperturba- 
tive long-distance contribution. In particular, not only the LCWFs for the light mesons n,p,..., 
produced in the final state, but also those for the B meson participate in processes where large mo- 
mentum is transferred to the soft spectator quark via hard gluon exchange ||l|]. We define the leading 
quark-antiquark component of B-meson LCWF as the vacuum-to-meson matrix element ^ : 

4>B{t,^) = (0|^(m)Pe'^^o'^^"-^(^"V75/2v(0)|B(v)) = I d(oe-''"'Mo},^)- (1) 

Here the bilocal operator is built of the Z^-quark and light antiquark fields, /Jv(0) and q{tn), linked 
by the Wilson line at a light-like separation tn with the light-like vector (n^ = 0, n-v = 1), is 
4- velocity of the B meson, and the RHS defines the momentum representation with G)v+ denoting 
the LC component of the Ught antiquark's momentum. A difference compared with the familiar 
pion-LCWF is that hy(0) is an effective field in the heavy-quark effective theory (HQET) [^]. In 
(|l|), n denotes the scale where the bilocal operator is renormalized. If one could Taylor expand the 
bilocal operator, the WF (|l|) would be equivalent to matrix element of a set of local operators of the 
type, q{0){n ■ D)-' ^Y5hy{0). Indeed, for the light mesons n,p,..., such correspondence holds and 
is useful for establishing systematic framework to calculate the LCWFs explicitly [Q]. In contrast, 
for the B-meson LCWF (|l]), such straightforward correspondence is broken as demonstrated below, 
because the twist is not a good quantum number by presence of the effective heavy-quark field. 

The "IR structure" of the B-meson LCWF has been studied using constraints from equa- 
tions of motion (EOM), qB = v- l^hy = 0, and heavy-quark symmetry, i^hy = hy. These constraints 
are solved to give (|l]) as the sum of the two terms, 0e(G)) = (l)jj^^\co) + ^jf\co). The first term 
is expressed by the analytic formula as ^^^^ (&)) = iFo)d{2A -co)/ (2A^), in terms of the HQET 
parameters [^], the decay constant F = —i{0\qijlY5K\B{v)) and the "effective mass" A = ms — mi, 
representing the mass difference between B-meson and Z>-quark. The second term 0^^^ is obtained 
as a certain integral of matrix element of 3-body LC operator, {0\q{tn)Gccp{un)a;t^r]Jh{^)\B{v)), 
where' the nonperturbative gluons generated inside the B meson participate as the field strength 
tensor G^jg. On the other hand, the "UV structure" of the B-meson WF has been studied in 
by calculating the 1-loop renormalization of the bilocal operator in ([T]). The "vertex-type" cor- 
rection around a "cusp" between the two Wilson lines, the light-like one of ([T]) and the time-like 
one from hy{0) = F exp[ig J^^dXv ■ A{Xv)]hy {—oov), produces the "radiation tail" as ^s(ft),/x) 
—iFCFOLsln{(o/^)/{K(o) for &)»/!, where Cf = (A'^ — l)/{2Nc)- This implies that the moments 
Jq dco(0^<pB{(0,IJ.), which would correspond to matrix element of the local operators mentioned 
above, are divergent reflecting the singularity associated with the cusp of the two Wilson lines 

To see the behavior of the B-meson WF incorporating both IR and UV structures, we calculate 
the radiative corrections taking into account hard and soft/coUinear loops. We calculate the 1-loop 
diagrams for the 2-point function (|l|) in 4 — 2e dimensions. We first note that the similar diagrams 
for the case of the pion LCWF, defined by the light-^^ bilocal operator, yield the "cancelling" UV 

' The Wilson lines connecting the constituent fields are suppressed for simplicity. 
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and IR poles, 1 /euv — l/sm, reflecting the scaleless loop-integral in massless QCD, and this struc- 
ture is accompanied by the convolution with the splitting function (Brodsky -Lepage kernel) [Q]; 
the corrections are analytic (Taylor expandable) at f = 0. Now, for the B-meson LCWF (|), we get 







(---) 
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+ •■•, (2) 



where L = ln(;YjLte^^) with the MS scale jj., and (•••) = (0| •• • \B{v)). The "vertex-type" correction 
that connects the light-like Wilson line and the q{tn) field in (|l]) is associated with only the massless 
degrees of freedom, and gives the term with the "plus"-distribution [(§/(! — ^)]+, which has the 
same structure as the corresponding correction in the pion case. But the other diagrams give novel 
behavior with "non-cancelling" UV and IR poles: the another vertex-type correction gives ||^ the 
terms proportional to 5(1 — which contain the double as well as single UV pole, corresponding 
to the cusp singularity mentioned above. The "ladder-type" correction, connecting the two quark 
fields in @, gives the last line in (|2|), which contains the IR poles and is associated with not only 
the bilocal operator in (|I|) but also the higher dimensional operators; note that the ellipses in (^ are 
expressed by those operators involving the two or more additional covariant derivatives. 

The renormalized LCWF, obtained by subtracting the UV poles from (^, is non-analytic at 
f = by presence of logarithms L, In particular, the nontrivial dependence on f/x through 

L implies that the scale ~ 1 plays a role to separate the UV and IR regions. Thus we have to use 
the operator product expansion (OPE) to treat the different UV and IR behaviors simultaneously: 
the coefficient functions absorb all the singular logarithms, while, for the local operators to absorb 
the IR poles, we have to take into account many higher dimensional operators. Such OPE with 
local operators is useful when the separation t is less than the typical distance scale of quantum 
fluctuation, i.e., when f < l//x. The OPE result can be evolved to a higher scale by the evolution 
operator associated with the single-UV-pole terms in (^), i.e., by the Brodsky-Lepage-type kernel 
and the Sudakov-type operator with the cusp anomalous dimension [^, ^. 

An OPE for the B-meson LCWF ([|) has been discussed in [^, taking into account the local 
operators of dimension less than 5 and the NLO {0{as)) corrections to the corresponding Wilson 
coefficients in a "cutoff scheme", where specifically an additional momentum cutoff Auv is intro- 
duced and the OPE is derived for the regularized moments, J^"^ d(0(0^ ^b{c^,1^) with j = 0, 1. In 
this work we derive the OPE for the B-meson LCWF (|I|), taking into account the local operators 
of dimension less than 6 and calculating the radiative corrections for the corresponding Wilson 
coefficients at NLO accuracy. Following the above discussion, we carry out the calculation in the 
MS scheme for f < l//x, so that there is no need to introduce any additional cutoff. 

The most complicated part of this calculation is to reorganize the contributions from a gauge- 
invariant set of Feynman diagrams in terms of matrix element of gauge-invariant operators includ- 
ing higher dimensional ones. In particular, to derive the Wilson coefficients associated with the 
dimension-5 operators such as qGapi^Jshy, we have to compute the 1-loop diagrams for the 3-point 
function, as well as those for the 2-point function as in (^, where the former diagrams are obtained 
by attaching the external gluon line to the latter diagrams in all possible ways. To perform the 
calculation in a systematic and economic way, we employ the background field method We de- 
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compose the quark and gluon fields into the "quantum" and "classical" parts, where the latter part 
represents the nonperturbative long-distance degrees of freedom inside B meson as a background 
field and satisfies the classical EOM. The quark and gluon propagators for the quantum part con- 
tain the coupling with an arbitrary number of background fields, and each building block of the 
Feynman diagrams obeys the exact transformation property under the gauge transformation for the 
background fields. We use the Fock-Schwinger gauge, x^A)i {x) = 0, for the background gluon 
field a[[^. This gauge condition is solved to give A^^\x) = /J duiix^ G^pj^{iix) [^, which allows us 
to reexpress each Feynman diagram in terms of matrix element of the operators associated with the 
field strength tensor. 

Using this framework, the tree-level matching to derive our OPE can be performed replacing 
each constituent field in (|I]) by the corresponding background fields. The classical bilocal oper- 
ator can be Taylor expanded, and we obtain the OPE at the tree level with the local operators of 
dimension-3, -4, and -5 (see (^) below). 

For the 1-loop matching, we calculate the 1-loop corrections to the 2- and 3-point functions 
with the insertion of the bilocal operator in (|l]), taking into account the mixing of operators of 
dimension less than 6. Apparently the mixing through the UV region of the loop momenta arises 
only in the 2-point function, whose result can be immediately obtained from @; but the additional 
mixing can arise in both 2- and 3-point functions accompanying the IR poles. We perform the 
loop integration in the coordinate space using the Schwinger ("heat kernel") representation of the 
Feynman amplitudes under the background fields [^. For the calculation of the 3-point function, 
the Fock-Schwinger gauge ensures that the Wilson line in (|T]), as well as the heavy-quark field, does 
not couple directly to the background gluons. We reorganize the result in terms of a complete set 
of gauge-invariant operators using the EOM and heavy-quark symmetry. We subtract the UV poles 
to renormalize the bilocal operator of (|l|), and also renormalize the local operators to absorb the 
IR poles. Combining this result with the above tree-level result, we obtain the OPE for the bilocal 
operator in q{tn)Pe'^^od^nMXn)^^^}^^(Q^ ^ ^,C;(f,/i)^,(;U), to the desired accuracy as 
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where Cg 

(3) ^ 
1 



Nc and a basis of local operators. 



as 



1,2,...), of dimension- J is defined 
q^Yshv, {^i^^} = {q{in ■ 'D)4ysh^, q{iv ■ ^)^75/j„}, and {fff^} = {q{in ■ 'Bfr^ysh,., 
q{iv-^){in-'D)4yshy, q{iv -"Df^y^h,., qigGa^v'^n^^y^h,,, qigOapfrfi^y^hy, qigGapfvP$y5h,,, 
^sGap'^'^^^Ys^v}- Here we have introduced another light-like vector, if' = 0, as = +n^)/2. 
The double logarithm in the coefficient functions originates from cusp singularity. The 1-loop 
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corrections for the 3-point function induce only the operators 



4,5,6,7 



associated with the field 



strength tensor, while those for the 2-point function induce all ten operators through the use of the 
EOM. 

Taking the matrix element (•••) = (0| •• • |5(v)) of @, we can derive the OPE form of the B- 
meson LCWF (|I]). Matrix elements of the local operators appearing in the RHS of (|3|) are known to 
be related to a few nonperturbative parameters in the HQET, using the EOM and heavy-quark sym- 
metry as demonstrated in [Q, ^]: = 4/FA/3, {^'^^) = iFA, where F and A were introduced 
below (|l|), and all seven matrix elements {^^{'^) for the dimension-5 operators can be expressed 
by f , A and two additional HQET parameters Xe and Xr, which are associated with the chromo- 
electric and chromomagnetic fields inside the B meson, respectively, as (^^E • ajshy.) 
{qgU ■ ajshy) = iFXji in the rest frame, v = (1,0). We get the OPE form for the LCWF (|l|) as 
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with the MS scale n and = cCs{lJ.), which takes into account the Wilson coefficients to 0{as) 
and a complete set of the local operators of dimension less than 6. Fourier transforming to the 
momentum representation and taking the first two regularized-moments, Jq"^ dO)0)^(j)B{(0,IJ.) (j = 
0, 1), the contributions from the first line in (Q), associated with matrix element of the dimension- 
3 and -4 operators, coincide completely with the result obtained in The second and third 
lines in are generated from the dimension-5 operators. Our OPE result @ "merges" the IR 
and UV structures peculiar to the B-meson LCWF ||5|, ^, and embodies novel behaviors that are 
completely different from those of the pion LCWF: /i and t are strongly correlated due to the 
logarithmic contributions from cusp singularity in the radiative corrections, so that the WF is not 
Taylor expandable about t = 0. The WF receives the contributions from many higher dimensional 
operators, in particular, from those associated with the long-distance gluon fields inside B-meson. 

Our result @ allows us to parameterize all nonperturbative contributions by a few HQET pa- 
rameters. As a result, the B-meson LCWF ([l]) obeys the two-step evolution: it is governed by the 
Sudakov-type and Brodsky-Lepage-type scale dependence from the high scale jXi ~ s/mb^^^QCD, 
associated with the QCD factorization formula, to the scale jj. < l/t, while that for the lower scale 
is governed by the anomalous dimensions of the local operators associated with the HQET pa- 
rameters. This result is useful for clarifying model-independent properties of the B-meson LCWF, 
combined with nonperturbative estimates of the relevant HQET parameters. 

K.T. thanks V. M. Braun for valuable discussions. This work is supported by the Grant-in-Aid 
for Scientific Research No. B-19340063. 
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